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Some recent progress in China in the study of charmless B decays with QCD factorization 
is reviewed. Chirally enhanced power corrections and infrared divergence problem are 
stressed. 

1. Introduction 

B meson weak decays are very important for testing the standard model and probing 
new physics. The two B factories have accumulated large data sets while Tevatron 
and LHCb will have even more data. The main task for theorists is to compute 
all experimental observables in a reliable way. This amounts to compute hadronic 
matrix elements reliably. There are several methods in the literature. The earliest 
is the so called navie factorization In this method the matrix element has no 
renormalization scheme and scale dependence, so cannot cancel the corresponding 
dependence of the Wilson coefficients. Futher more, it cannot account for "non- 
factorizable" contributions. The generalized factorization El takes into account the 
one-loop radiative corrections to recover the renormalization scheme and scale de- 
pendence for the hadronic matrix elements. But at quark level, to avoid the infrared 
divergence, it has to assume the external quarks to be off-shell by — p 2 . This results 
in gauge dependence. To account for "nonfactorizable" contributions, a phenomeno- 
logical parameter N^ s is introduced and it is assumed that N° s is universal. But 
actually N^ s is process-dependent. 

In order to overcome the shortcomings of the above methods, in 1999, Beneke, et 
al. (BBNS) E l proposed a new scheme based on QCD. In this talk, I will concentrate 
on our work rM5 | 6 | 7 | c >J aDOU t the chirally enhanced power corrections and infrared 
divergence problem and the application of the QCD factorization approach. My 
talk is organized as follows: In Section 2, I shall give an simple introduction to 
QCD factorization. Section 3 is devoted to the chirally enhanced power corrections 
and the cancellation of the infrared divergence in vertex corrections. In section 4^6, 
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the application of QCD factorization including chirally enhanced power corrections 
to the two-body charmless £?-decays is reviewed. Section 7 is for conclusions. 

2. General remarks on QCD factorization 

Beneke, et al. (BBNS) proposed a new factorization scheme based on QCD They 
pointed out that, in the heavy quark limit (m^ — > oo) 

pi 

(MiM 2 \Qi\B) = Ff^ M2 (q 2 ) dxTf(x)* Ml (x) + (Mi~M3) 

J o 

+ V / didxdyTl J I {i,x,y)^ B {^M l {x)^ M M + O^QCD/m b ) (1) 
. J 

where the short distance hard scattering kernels T 1 ' 11 is calculable order by or- 
der in perturbative theory; the long distance quantities, e.g. decay constants, form 
factors, light-cone distribution amplitudes are inputs form either experimental mea- 
surements or other theory; M\ is a light meson or a charmonium state, Mi (contains 
the spectator in B) is any light or heavy meson. If Mi is heavy (for example D), the 
second line in Eq.lQJ is l/m& suppressed. At the zeroth order of a s , Eq.Q reduces 




(a) (b) (c) (d) 




(e) (f) (g) (h) 

Fig. 1. T 1 ' 11 at 0(a s ). The upward quark lines denote the ejected meson Mi 



to naive factorization. At the higer order, the corrections can be computed system- 
atically. The renormalization scheme and scale dependence of {Qi) is restored. In 
the heavy quark limit, the "nonfactorization" contributions is calculable perturba- 
tively. It does not need to introduce N° s ■ The strong phase is suppressed by either 
a s or l/mft. W^-exchange and VF-annihilation diagrams are l/m& suppressed. There 
is no long distance interactions between Mi and (-BM2). T/ includes tree diagram, 
non-factorizable gluon exchange, e.g. vertex corrections [Fig^a)-(d)], penguin cor- 
rections [Fignje)-(f)], ' ' '• Tfj includes hard spectator scattering [Fig^g)-(h)], • • • 
Later on we shall apply Eq.QJ to B meson rare decays. 



2, 2008 3:51 WSPC/INSTRUCTION FILE august 



QCD Factorization and Rare B Meson Decays 3 

3. Chirally Enhanced Power Corrections and Infrared Divergence 
Cancellation in Vertex Corrections 

Superficially, ^QCDjfnb ~ 1/15 is a small number. But in some cases, such power 
suppression fails numerically. For example 

(Q 6 ) F = -2j2{Pi\(m')s-p\0){P2\ms +P \B) (2) 
i 1 

is always multiplied by a formally power suppressed but chirally enhanced factor r x 
= 2jj,p/mb (where \ip = m P /(m\ + 7712), mi,2 are current quark masses, and P1.2 
denote pseudoscalar mesons). 

In the heavy quark limit, rrib — > 00, r x — > 0. But for ~ 5GeV, r x ~ C(l), no 
l/m^ suppression. So we should consider the chirally enhanced power corrections. 
For example, in B — > 7r.ff, dominant contribution to the amplitude ~ + a6r x , 
and 04 ~ a6 r x' so one na ^ comes from the chirally enhanced power corrections. 

Possible sources of power corrections are : high twist wave functions, quark 
transverse momentum k±, and annihilation topology diagrams. In this talk, I only 
discuss the chirally enhanced power corrections from high twist wave functions ^1 

All the chirally enhanced power corrections involve twist-3 light cone distribu- 
tion amplitudes 3> p (x) and If we want to include chirally enhanced power 
corrections consistently, we must prove that the hard scattering kernels are infrared 
finite, at least for vertex corrections. 

The decay amplitude of B — ► M\Mi reads 

A(B^M 1 M 2 ) = ^= E^ a i(M)(^iM 2 |Qi|5) F (3) 

V q={u,c} i=l 

For the hard kernel T? in Eq. Q , there is no contributions from twist-3 wave function 
insertion for (V - A)(8>(V - A) and (S + P)(8>(£ - P) operators. But for (V + 
A)®(V — A), it is subtle! For twist-3 light cone distribution amplitudes, if we use 
dimensional regularization, the infrared divergences of Fig.l(a)-(d) can not cancel! 
This is because the twist-3 wave functions arc only defined in 4 = 3 + 1 dimensions. 
We need to give gluon a small mass and regularize the infrared integrals in 4- 
dimensions. 

After a lengthy calculation, for the operator Q5, the vertex correction terms are 



a s C F § a {v) r ln 2 A 
4tt N c v i 2 

Fig. 1(c) ~ + — £lL^M. hn 2 a - 21n(-«)lnA + 31nA + finite terms! (6) 

4-7T Nc V 1 J 

Fig.l(d) as_CF^v) T ln 2 A _ 2 \ n (-v)hi\ + 31nA + finite terms! (7) 

4-7T JV C v L J 



Fig.l(a) — — \ l JL2l + 21n(— u)lnA - 41nulnA + InA + finite terms! (4) 

47r JVc v L 2 J 

Fig.l(b) ~ + ^l^f^M. \ l JL* + 21n(-u)lnA - 41nt4nA + InA + finite terms! (5) 



where A = m^/m^. 
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In the above, we can see that after summing over Fig^a)-(d), the infrared 
divergences cancel only when ^ a {v) = <& a (v), where V = 1 — v. This is not a surprise 
because we neglected the contribution of three-body twist-3 light cone distribution 
amplitudes (LCD As). For the asymptotic form, 3? CT (i>) = 6vv which is the same for 
both 7r and K meson. 

After our work^l BBNS ^ also discussed the infrared divergence problem. They 
implemented the equation of motion 

fi(fci) jk 1 ---v(k 2 ) = (8) 

to show the infrared safety without assuming 3> ff (u) to be symmetric. However 
Eq.(JH| is justified only when 2-particle twist-3 LCD As are asymptotic, i.e. & p (v) 
— 1 and <fr a (v) — 6vv by neglecting the contribution from the 3-particle twist-3 
LCDAs. Our result is based on a more general situation than that BBNS considered. 
Therefore, BBNS result is consistent with ours. 

We can prove at the order of 0(a s ) that the decay amplitude is independent 
of the renormalization scale. We can also prove the gauge invariance when chirally 
enhanced power corrections are included. The detailed proof can be found in El We 
do not discuss it here. Now we come to the application below. 



4. B — > PP, PV Charmless Decays and CP Violation 

Now we give the prediction of the branching ratios and CP asymmetries. We in- 
cluded the contributions of the hard spectator scattering and annihilation topology. 
We also include the chirally enhanced power corrections. We compute the branching 
ratios (Br) and CP asymmetries of B — > PP : tttt, nK, KK, Kr)('\ irrf'\ rf'^rf'^, ■ ■ ■ 

■ ■ ■■ B -> PV : up, ttuj, nK* pK , cuK, K*^'\ pjjW, .All our numerical results 

can be found in the tables in Refernce'^l From our calculated resulte we see that : 
Owing to including the chirally enhanced power corrections, the scale p dependence 
of Br is smaller; Br of S ± — > tt° tt~ (pure tree) and B ± — ► X°7r ± (pure penguin) 
are in good agreement with data; Br of B° — > ir + ir~ is larger than data; Br of B 
— > nK seem smaller than data; For B — > Kri^\ if consider di-gloun fussion, we can 
fit the data; For B° — > 7r°7r°, our prediction is much smaller than data; For B° — > 
K + K~ , the result is very small, only weak annihilation diagram contributes; There 
are large uncertainties from CKM elements, form-factors, annihilation parameters. 
But a global analysis can fit the data very well. 

For CP asymmetries, the numberial results are not reliable because vertex, pen- 
guin, hard spectator scattering and annihilation diagrams can all give imaginary 
part to decay amplitudes, so strongly affect Acp- Because of space limitation, for 
more detail, see 1^1 For B — > PV, the differences from B — > PP arc: If the emitted 
meson is vector meson (vertex diagrams and penguins), then the twist-3 wave func- 
tions of V do not contribute (power suppressed, so can be neglected) . But for hard 
spectator scattering, even the emitted meson is vector one, there are still twist-3 
contributions. Note, there are large uncertainties on CP asymmetries. So we cannot 
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make good predictions on Acp in QCDF scheme. The calculated branching ratios 
of B —> PV for b — > d and b — > s transitions and CP asymmetries can be found in 
Reference "4 



5. Global Analysis of B — > PP, PV Charmless Decays 

Beneke, et al. ™ have done a global analysis on B — > PP. But their fitted 7 ~ 90°, 
a bit large compare with the standard CKM fit. Now there are many new data on 
both B -> PP and PV. It is necessary to do the global fit of B -> PP, P^ at the 
same time. We have done it! 

For P -> tttt, irK, the fit is sensitive to \V ub \, j{p,r)), F B ^ n , P B ^ K , X A , 
}b/^b, and m s . To include seven P — > PV decay channels, only Aq^ p , X BV are 
newly involved sensitive parameters. So the inclusion of B — ► PV^ will lead more 
stringent test to QCDF. The experimental constrains of CP asymmetries are not 
implemented, because the QCDF's predictions on CP asymmetries in B decays are 
rough! So we do not use CP asymmetries for the global fit. The best fit values of 
the branching fractions is presented in Tabled 



Table 1. Fitl and Fit2 mean the best fit values of Br in unit of 10 — 6 with and with- 
out the contribtuions of the chirally enhanced hard spectator and annihilation topology, 
respectively. 





modes 


Exp. 


Fitl 


Fit2 




modes 


Exp. 


Fitl 


Fit2 




— > 


4.77±0.54 


4.82 


5.68 


B+ 


tt+tt" 


5.78±0.95 


5.35 


3.25 


B° 


-> K+ir- 


18.5±1.0 


19.0 


18.8 


B+ 


-> K+n° 


12.7±1.2 


11.4 


12.6 


B+ 


-> K°TT + 


18.1±1.7 


20.1 


20.2 


B° 


-> X°7r° 


10.2±1.5 


8.2 


7.3 


B+ 


— ► JJ7T + 


< 5.2 


2.8 


1.8 


B° 


-» 7r±pT 


25.4±4.3 


26.7 


29.5 


B+ 


-> 7T+p° 


8.6±2.0 


8.9 


8.5 


B° 


-> if+p" 


13.1±4.7 


12.1 


5.1 


B+ 


-* 4>K+ 


8.9±1.0 


8.9 


7.1 


B° 


-» </>K° 


8.6±1.3 


8.4 


6.7 


B+ 


-> r;p+ 


< 6.2 


4.6 


3.8 


B° 


-» o,K° 


5.9±1.9 


6.3 


1.2 



The main results are : 

The best fit of 7 ~ 79° which is consistent with recent fit results 37° < 7 < 80^. 
For B -> 7r°7r°, the best fit is around lxl0~ 6 , 

Exp: the BaBar and Belle average is (1.90±0.49)xl0 -6 
For B+ -> u)K + , the best fit is 6.25xl0~ 6 , 

Exp: Belle (9.2t^±1.0) x 10~ 6 CLEO < 8xl0~ 6 . 
For B+ -> LO-K+, the best fit is 6.66xl0~ 6 , 

Exp: BaBar (6.6l^±0.7)xl0" 6 Belle < 8.2xl0~ 6 
For B+ -> tt+K*°, the best fit - 10xl0~ 6 , 

Exp: LP03 - (10.3±1.2t^)xl0- 6 ICHEP04 - (9.0±1.3)xl0" 6 . 
For more details, see Reference^. 
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6. B s —> PP, PV Charmless Decays 

We include : i) chirally enhanced power corrections, ii) weak annihilation, iii) hard 
spectator scattering. We list the computed branching ratios and CP asymmetries of 
B s -> PP, PV in Referenced Only B s -> K&K (HT 6 ~Hr 5 ), K^if* (10~ 6 ), 
K ± p T (10~ 5 ), rf'^rfi'} (10~ 6 ~10~ 5 ) have large branching ratios. 

7. Conclusions 

We have shown that 

• chirally enhanced power corrections must be included in QCDF 

• twist-3 light cone distribution amplitude <& a {x) must be considered 
simultaneously 

• the infrared divergences in vertex corrections cancel only twist-3 wave func- 
tion of light pscudoscalar is symmetric, so chirally enhanced power en- 
hanced power corrections can be included consistently 

• the calculated branching ratios for B — > PP, PV charmless decays are, 
principally, in agreement with data 

• gloabl analysis of B — > PP, PV can fit the data very well including B° — > 
7r°7r° and 7 - 79° 

• B s — > PP, PV charmless decays are waiting to be tested in Tevatron and 
LHCb 

I thank H. J. Gong, J. F. Sun, D. S. Yang, and G. H. Zhu for their collaborations 
and discussions. 
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